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1. Introduction 

High energy QCD has reached a mature stage of development, for dilute-dense scattering (for example for DIS 
with nuclei) [1, 2, 3, 4, 5, 6]. However, for dilute-dilute scattering at high energy, (for example the scattering of 
two virtual photons with large but almost equal virtualities), despite the great deal of effort by experts in the field 
(see for example Refs. [7, 8, 9, 10, 11, 12, 13, 14, 15] ), we are still waiting for the desired breakthrough. In this 
paper we address the problem of the scattering amplitude for the dilute - dilute system, using the BFKL Pomeron 



calculus [16, 30, 1, 2, 12, 16, 17]. The BFKL Pomeron calculus is elegantly formulated in terms of the functional 
integral [12]. Namely, 

Z[$,$+] = J D^D^+e^ with S = Sq + Si + Se (1.1) 

where Sq describes free Pomerons, Si corresponds to their mutual interaction while Se relates to the interaction 
with the external sources (target and projectile). Sq has a simple expression: 

So = J dY dY' (fxi(fx2(fx[(fx2<P^{xi,X2;Y)G-^{xi,X2;Y\x[,X2;Y')<^{x[,x'2;Y') (1.2) 

while 



Si='^ [ dY' [ ^!^^?p^|(Li2$(xi,X2;y')) ^-'{xi,xs;Y')^+{xs,X2;Y') + h.c.} (1.3) 

iVc J J X'^2 ^23 -^IS 

where h.c. denotes the Hermitian conjugate. Se describes the interaction with the scattering particles (two 
small colliding dipoles in our case), but we do not need the explicit expression for this term in our paper. Y is the 
rapidity of the dipole. For dipolc-dipole scattering 1" = In s in the leading log approximation where s is the energy 
of colliding dipoles. Eq. (1.1) is written in the leading log approximation. The following notation, = Ncas/ir will 
be assumed throughout. $ and relate to the BFKL Pomeron and G{xi,X2;Y\x[,X2;Y') is the Green function 
of the BFKL Pomeron which takes the form: 

G-\xi,X2;Y\x[,x'2;Y') 

with 

Hf{xi,X2;Y) = ^ J d'^xsK {xi,X2\x3) {f{xi,X2;Y) - f{xi,X3;Y) - f{x3,X2;Y)} (1.5) 
where 
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K{xi,X2\x3) = 2^'^2 -^12 = x\2piP2 Pk = -i^ xu = 1> 2) (1.6) 



The exact Green function for the Pomeron, which is the goal of this paper to derive, is equal to: 

G{x^,X2;Y\x,,X2;Y) = / e^[^.^+] ^^'^^ 



and the bare (initial) BFKL Pomeron Green function [30] is determined by Eq. (1.7), where only the term S'o is 
included in S [$,$+]. 

It is worthwhile mentioning that due to conformal invariance of the BFKL Pomeron calculus, the form of the 
triple Pomeron vertex is known unambiguously, and coincides with the direct calculations found in ref. [17]* The 
simplicity of Eq. (1.1) is thanks to the key assumption, that the triple Pomeron vortex is the only essential vertex 
for Pomeron interactions. In other words, in Eq. (1.1) we neglect all other local vertices (for example, the vertex 
for the transition of one Pomeron to three and so on). We have no rigorous proof for this conjecture. There even 
exist arguments that the four Pomeron vertex should also be included [18]. It should noted also, that the BFKL 
Pomeron calculus given by Eq. (1.1), is formulated in the leading \/Nc approximation, where Nc is the number 
of colours. It is known that the BFKL Pomeron calculus cannot be a correct approximation for the scattering 
amplitude in the next-to-leading order in 1 /N^ approach, due to the fact that 2n- gluon states in the t-channel give 
a larger intercept than n BFKL Pomerons [19]. Nevertheless we consider Eq. (1.1) to be a good first approximation, 
for the dilute-dilute system of scattering. We will return to the discussion of all these problems in the conclusion, 
where we show that they are not important for the solution derived in this paper. 

The goal of this paper is to calculate the class of enhanced diagrams. The simplest examples of these diagrams 
are shown in Fig. 2 and Fig. 5. In other words we are going to sum BFKL Pomeron loops in this paper. These 
diagrams lead to a new Green function of the BFKL Pomeron (the term in Eq. (1.1) ), while the vertices of the 
interaction of the new dressed Pomeron remain the same as they appear in Si . It is worthwhile mentioning that if 
the Green function of the dressed Pomeron, will be such that the Pomeron contribution will lead to a decrease with 
energy, then the problem would be solved without needing to consider the interaction of the dressed Pomerons. 
However, if this is not the case and the dressed Pomeron still increases with energy, then the interaction of the 
dressed Pomeron needs to be included. On the other hand, for the scattering of two dipoles with small sizes, within 
a wide range of energy the enhanced diagrams dominate, since the interactions of the Pomeron with the target and 
the projectile are small in the leading 1/iVc-approximation. 

This paper is organized in the following way. In the next section the formulae for the triple Pomeron vertex 
and the Pomeron Green function are introduced. Using these ingredients we calculate the sum over the class of 
enhanced diagrams shown in Fig. 2 - Fig. 5, namely one-Pomeron loops in scries, in (cj, u) representation. This is 
done in a step-by-step way, in order to introduce the reader to the method of integration and assumptions used in 
this treatment, which will be used later on for more complicated diagrams. 

The third section is the main body of this paper. Using the techniques developed in section 2, we extend 
this approach to the sum over all enhanced diagrams, using two principle selection rules. First, we are searching 
only for the contribution to the vertices that are singular in v. Second, we assume that the contribution of the 
Pomeron loops in effective vertices, for multi-Pomeron production are negligibly small. This approximation is 
closely related to the Mueller-Patel-Salam-Iancu approach [20], but in this paper this strategy is formulated in 
(a;, z/)-representation. In section 3.1 the equations for the effective multi-Pomeron vertices are derived and solved. 

*For the sake of completeness in this presentation, we would like to mention that there is still a discussion in the literature, about 
whether or not the BFKL Pomeron calculus in the form of Eq. (1.1), correctly takes into account the reggeized gluons. However, Eq. 
(1.1) reproduces the non-linear Balitsky-Kovchegov equation, including the term for the gluon reggeization. As far as we know, no other 
examples have been suggested, where Eq.(l.l) provides an incorrect result due to problems with gluon reggeization." 



Figure 1: The basic diagram with the exchange of one Figure 2: The diagram with one Pomeron loop. 

Pomeron in the t-channel. 



In section 3.2 we present the results of the summation, for the Green function of the dressed Pomeron. It turns 
out that the exchange of the dressed Pomeron, leads to the total cross section for dipole-dipole scattering that 
decreases with energy. This result is reminiscent of the Pomeron Green function in 1 + 1 dimensional Pomeron 
calculus, [21, 22, 23, 24, 25] which also decreases due to the contribution of Pomeron loops. In the conclusion 
section, we summarize our results, and discuss the current stage of development of the BFKL Pomeron calculus. 
The appendix provides detailed information about the ingredients of the BFKL Pomeron calculus, including the 
bare Pomeron Green function and the triple Pomeron vertex. 

2. Summing simple Pomeron loops 

In this section we sum the set of enhanced diagrams shown in Fig. 2 - Fig. 5, namely one-Pomeron loops in series. 
The sum over all diagrams of this type will teach us the main characteristic features of the enhanced diagrams, 
which we will use for developing an approach for summing over a more general class of diagrams. 

2.1 Bare Pomeron 

The most basic diagram is the exchange of one Pomeron in the t-channel without any loops. The expression for 
this diagram which is shown in Fig. 1 is given by the following expression in oj representation [26, 27, 35, 28, 29]: 



,2 



UJ — UJ {v) 



where g {v) = 




2 



(2.1) 



The following inverse Mellin transform allows us to pass to Y representation: 




(2.2) 



The contour of integration is placed to the right of all singularities of Gq {ujv). 



2.2 One loop 



The first correction to tlic basic diagram is the diagram witli one loop shown below in Fig. 2. The amplitude for 
the one-loop diagram takes the following form in a; representation [12]: 



A(^i^{ui,i') = Go(w,z^) / ciz^' m (w, z^, I/') Go (w, z^') (2-3) 

J — oo 

where m{i',i/',u) is the Pomeron self mass, given by the following formula: 



/•oo roo r I U 

m{u,u',u) = dui du2 g {vi) g {v2) ■ 



-oo •/ — oo 



U — UJ {Vi) — U {U2) 



(2.4) 



The triple Pomeron vertex has the following definition [34] ^: 

r {v I VI, V2) = (^rpianar {'^ I ^l, ^2) - ^non-p\&nw I ^2)^ (2.5) 



rplanar(2^kl>'^2) = (^4 + ^ J Q.{y\vi,V2) (2.6) 

rnonplanar(i^kl,i^2) = (4+^^) A (^^ ^1, 1^2) (^X (i^) - X (l^l) - X (1^2)^ (2.7) 

where the functions X ('^) ) ^ (^^ | ^^i, 2^2) and A (z^ | z^i, 1/2) are defined explicitly in Eqs. (A. 2), (A. 6) and (A. 24) 
respectively. rpia,nar and ^riorLp\&nax denote the two diagrams that contribute to the vertex, namely the planar and 
non-planar diagrams shown below in Fig. 3 (a) and Fig. 3 (b) (the diagrams in Fig. 3 are taken from ref. [34]). 

r {y I z^i, ^2) is the vertex at the top of the loop of Fig. 2, whereby the Pomeron with scaling dimension u splits 
into two Pomerons that form the loop, with scaling dimensions v\ and V2- The vertex at the bottom of the loop is 
labeled by F (z^' | i^i, zv2)) whereby the Pomerons with scaling dimensions vi and 1^2 recombine to one Pomeron with 
scaling dimension u' . Before proceeding to calculate m (w, z/, z/'), it is instructive to switch to the variables A, a 
and A defined as: 



\ = i/2-v; A' = z/2-z/'; cr = z^i + 1'2 ; A = z^i - 1/2 (2.8) 

^The coefficient in front of Eq. (2.5) contains an extra factor of 1/iVc in accordance Eq. (1.3), compared to the same coefficient that 
appears in [34]. 
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1 '2 
(a) The planar diagram contribution to the triple Pomeron vertex. 





1 2 
(b) The non - planar diagram contribution to the triple Pomeron vertex. 

Figure 3: The two contributing diagrams to the triple Pomeron vertex. 



In terns of (A, cr, A) variables, then Eq. (2.4) reads: 



m 



J —oo J —oo 



dAg 



a + A\ /a-A 



r X 



a, Ajrf A' 



a, A 



u — UJ 



a + A 



a-A 



(2.9) 



The product of vertices F (A | a, A) F (A' | a, A) that appears in Eq. (2.9), contains tv^^o poles in the a-plane in 
the following regions (for details see Eqs. (A. 16 - A. 18) and Eq. (A. 24)): 



a = X 



a = i — X 4^ 1/2 + 1/ = U1 + U2 



(2.10) 
(2.11) 



We need to check that the vertices and A converge as — >■ 00, for both of the regions of Eq. (2.10) and 
Eq. (2.11). Prom Eq. (A. 24) its clear that A falls down for pure imaginary u = ik, as exp (— 2Kln«;). Actually, 



we will see below that we need to integrate the z/-image of the amplitude to calculate the scattering amplitude in 
coordinate space, in the following way: 



^(0 = / ^e-"^<^H (2.12) 

One can see, that independent from the sign of In^, we have to close the //-contour of integration for A over 
the upper half-plane. Since the pole of Eq. (2.10) lies in the upper half-plane, then the //-integral for the part of 
the vertex proportional to A, is taken using the contour which is closed around the pole given by Eq. (2.10). 



The situation with the part of the vertex proportional to Q is quite different. One can see that at large and 
pure imaginary u, then U falls down only as a power of u. Therefore, the convergence of the integral in Eq. (2.12) 
depends on the sign of InC- We choose C > 1 • For this choice of ( we can close the //-integration contour over the 
upper half-plane for both A and il. 



With this in mind, A ^ from above the real axis, is our region of interest for calculating the contribution 
to the vertex proportional to A. Whereas A ^ from below the real axis, is the relevant region for calculating 
the part of the vertex proportional to fl. The dominant high energy behaviour exp (2a;o Y) stems from the region 
I'l, 1^2 — >■ 0, (or in other words a 0). Thus since A — )■ 0"^ is our region of interest, then Eq. (2.10) is the relevant 
region, whereas Eq. (2.11) does not yield a ^ as A — )■ 0"*". In light of this observation, it is useful to extract the 
pole explicitly by making the following definition: 



r(A|.,A) = mUl^ (2.13) 

A — (7 

In this approach, T (A | a, A) is finite in the region a — > A. Hence plugging Eq. (2.13) into Eq. (2.9) leads to 
the following formula: 



m (cj. A, A') 



^oo poo 

da / dA 

-oo J — oo 



V + A\ /a-A 
9\ ^^]9 



r A 



C7,A f A' 



cr, A 



UJ — OJ 



cr A 



a- A 



A-C7 A'-a 



(2.14) 



J — c 



-2m I dA ( 

-oo 



'a + A\ ( a - A 



r A 



c7,Ajf ( A'l cj,A 



LiJ — bJ 



a ^ A 



a- A 



A'-A 



(2.15) 



<T = A 



where in the last step the cr-integral was solved by taking the residue of the pole at o" = A. Switching back 
to (//, i/i,//2) variables, then Eq. (2.15) reads: 



m 



/oo 
dl'2 < 
oo 



g{vi)g{i^2) 



Ui,U2 f ( I/' 



1^1, 1^2 



U — UJ {vi) — OJ (^2) 



V — V 



(2.16) 



i/2—i/ = 1/1+1/2 



Interestingly, the Pomcron self-mass m (lv, v^v') contains a simple (first order) pole in the region v = v' . Hence 
using the result of Eq. (2.16): 



dv'm (cj, z/, z/') Go (w, z^') = 47r^ / dv^ < 



UJ — LO (l>\) — OJ (^2) 



Go (w, z^) 



(2.17) 



i/2— i/ = z/i+i/2 



where the right hand side of Eq. (2.17) was derived by taking the residue of the pole at = f' , after integrating 
over u' . The right hand side of Eq. (2.17) can be re-written in the following equivalent form: 



/oo 
dv' 
00 



m (a;, z/, z/') Go (o;, i^') = E (w, z^) Go (w, z/) 



(2.18) 



S (w, v) 



dvi / 

00 J —00 



dvo 



g [vij g [1^2) 



Z/|z/i,Z/2 



UJ — CO (Vi) — UJ (z/2) 



u — V1 — V2 



(2.19) 



where the vi integral in Eq. (2.19) is solved by taking the residue of the pole at i/2 — z/" — z/"! — z/'2 = 0. Thanks 
to the simplification of Eq. (2.18), then Eq. (2.3) can be re-cast as: 



^(1) (w, v) = Go (w, v) S (w, v) Go (w, v) 
In order to pass to Y representation, use the following inverse Mellin transform: 



(2.20) 



^(1) {Y,v) 



2m 



'o+«oo 



du e'^^74(i) (uj, u) 



(2.21) 



which after inserting Eq. (2.20) yields: 



Figure 4: The diagram with two Pomeron loops. 



Figure 5: The diagram with n Pomeron loops in suc- 
cession. 



g [1^1)9 [1^2 



= / di^i I du2- ^-^^ni^ -TT^ -\- (2.22) 



-oo 



-oo I a; — w (z^i) — oj (1/2) J { 2 ~ ^ ~ ^1 ^ ^2 



+ Y 



CO (u) — U! (vi) — U (^2) \ CO (ly) — 00 (vi) — CO {1^2) 
2.3 Two loops 

The second correction to the basic diagram is the diagram with two loops shown below in Fig. 4. The amplitude 
for the 2-loop diagram takes the following form in lu representation: 



A(2) (w, z^) 



/oo /*oo 
du' j du" m{ijj,u,u'^ Gq{oj,v') m (^uj,v' Go^uj,!/") (2.23) 
00 J —00 

Thanks to the useful result of Eq. (2.18), the amplitude for the 2-loop diagram of Eq. (2.23) simplifies to the 
following formula: 

^(2) (to, v) = Go {uj, v) ( Go {to, v) S (w, v) \ (2.24) 



where E(a;,z^) is given by Eq. (2.19). In order to pass to Y representation, use the following inverse Mellin 
transform: 



ra+ioo 

A(2)(F,i/) = ^ da; e'^^^(2)(u;,z^) (2.25) 



2m 

^ a—ioo 



"oo 



9'{u)ll / duU du^ g[,y^) g[4)j—^ ^ (2-26) 



fc=l — oo J — oo ( — — I' — I'l — 1^2 

wi Y UJ2 Y 

+ - + e-M^ 



I a; — a;i 00 — 102 \ ('^ ~ ^1) ~ ^2) / 2! 



du \{oj - uji) {uj - UJ2) J duj"^ \{uj - Ui) {oj - 0J2) J "^^^ 



ojk = oj(y\)+oj[4) (2.27) 

2.4 n loops 

Extending this approach to the diagram with n loops in succession shown in Fig. 5, leads to the following amplitude 
in oj representation: 

A(r^) {oJ, u) = Go {oj, u) Gq {uj, u) S {iO, v) ^ (2.28) 

The sum over the class of diagrams shown in Fig. 5 with an alternating minus sign, for all n G (0, oo) i.e. from 
n = loops up to infinity gives the Green function of the dressed Pomeron, labeled G2 [y-iV)- In this notation: 



G2(a;,^) = f;(-l)"A(„)(a;,^.) = ^^^^ ^ (2.29) 

n=o ( ^ + ^0 ^) ^ ^) 

From Eq. (2.29), the renormalized propagator G2 (w, v) can be expressed in terms of the bare Pomeron propa- 
gator Go (w, v) as: 

G-i(a;,i/) = G^^{uj,v) + S(a;,z/) (2.30) 



2.5 Calculation of S (a;, u) 

In this section, we calculate explicitly the formula for the Pomeron loop E (a;, u), defined above in Eq. (2.19)^ . In 
(A, (7, A) notation introduced in Eqs. (2.8): 




(2.31) 



(2.32) 



where in the last step, the integration over a was solved by taking the residue of the pole at cr = A. On 
the RHS of Eq. (2.29), there are singularities at i/2 + i/ — >■ and i/2 — u ^ that stem from Go (w,!/) in the 
numerator (see definition of Eq. (2.1)). However as discussed above, only the region i/2 — v = A — )■ is relevant 
(see Eq. (2.10) and the surrounding discussion). In light of this, the largest contribution to the propagator of the 
dressed Pomeron G2 {td, v) stems from the region i/2 — v = A — >■ 0. Hence we need to know the asymptote of 
E (w, v) in this region which can be found from Eq. (2.32): 



A + A\ /A-A 



lim S(a;,A) = -2m lim / dA ,y ^ - - f2 A|c7,A (2.33) 



Throughout the paper the notation \\m.x-^af{x) is used to denote the asymptotic behaviour of f(x) at a; — >■ a. It 
is not meant in the conventional way as used in analysis It is worthwhile mentioning here the two contributions to 
the vertex P (A | cr. A) in this region. Recall that from the definition of Eq. (2.5) and Eq. (2.13) that: 



f(A|a,A) = ^(A-a)(^Pplanar(Ak,A)--^r,o,planar(A|a,A)) (2.34) 

Prom Eqs. (A. 22) and (A. 26) the asymptotes (A — a) Ppianar (A | A) and (A — a) Pnonplanar (A | A) in the 
narrow region where A — >■ o" and o" — > 0, are the same up to a numerical coefficient. With this in mind, thanks to 
the suppression factor of 2'k/NI in front of Pnonplanar (A | cr. A) in Eq. (2.34), the contribution of the non-planar 
diagram to the vertex is parametrically smaller by a factor of 2tt/N'^ than the planar diagram, in our region of 

^Throughout the paper the notation linix-^a f{x) is used to denote the asymptotic behaviour of f(x) at a; — > o. It is not meant in the 
conventional way as used in analysis 



interest. Nevertheless for the sake of completeness, the contribution of both diagrams to the vertex are included in 
the calculation of S (w, A). Now inserting the asymptote of Eq. (A. 29) into Eq. (2.33) leads to the following result: 



lim S(w,A) 



-27ria^ lim / dA 
A^O 



A + A\ /A-A 



A + A\ /A-A 

[ U — iO { — UO 



\ 



A+A A-A 



(2.35) 



where the numerical coefficient a is given in Eq. (A. 28). Assuming that the typical value of A is small, then§ 
Eq. (2.35) can be re-cast as follows: 



lim T,(uj,X) = -27rza^lim / dA 

A^O A^O 



X + A\ /A-A 
9 



( 



-oo oj' 



A - A . A - A_ 



A+A A- A 



(2.36) 



(2.37) 



After closing the A integration contour over the upper half -plane that encloses the pole at A = A_|_, and 
taking the residue in this region, then Eq. (2.36) simplifies to: 



lim E(a;, A) 
A^O 



= AiT^a^ lim 



^ 'A + A+A /A-A, 



A. - A. 



A + A+ A - A 



= 27r^a^ lim 
A-)-0 



A + A, 



A- A_ 



uj -2uj [ - 



A + A+ A - A 



(2.38) 



Assuming that A+ is small (i.e. in the region a; — >■ 2a; (A/2)), then as A ^ 0: 



^Throughout the paper the notation linix-^a f{x) is used to denote the asymptotic behaviour of f(x) at x —¥ a. It is not meant in the 
conventional way as used in analysis 



' A + A 

lim g{ -^)^4(A + A+)^ (A+<1) (2.39) 



A^O V 2 

where the definition of Eq. (2.1) was used, with a similar result for ^ ( (A — A+) /2 ). Hence in the narrow 
region that oj ^ 2oj (A/2) and A ^ 0, then Eq. (2.38) reduces to: 



\4 

lim E(a;,A) = 327r2a2 lim = (2.40) 

bJ^2oj{\/2) w-)-2w(A/2) Joj - 2a; ( - 



for small A+, as a; — >■ 2a; (A/2). 



2.6 Green function of the dressed Pomeron 



The Green function of the dressed Pomeron can be calculated using Eq. (2.30), which can be reduced to the 
following expression 



G2^{u,v) = ^('^ - '^{y) + 5 (i^) S (a;, z/)) - [u - u {y) - (l/4A^) S (a;. A)) (2.41) 

The singularities of the Green function stems from the following equation (substituting Eq. (2.40)): 'plitude. 
The second is that they contribute in the rather narrow region A ^ and uj — >■ 2a;(A/2). 



'Throughout the paper the notation linix-^a f{x) is used to denote the asymptotic behaviour of f(x) at a; — > o. It is not meant in the 
conventional way as used in analysis. 

w — w (z/) — (l/4A^) E(w,A) = (2.42) 
a;^2w(A/2) - 2w 



u - - lim = = (2.43) 



2 

One can see that in the region where a; > 2w (-1), the correction to the pole at w = w(z/) is small, and can be neglected. However 
when a; — >^ 2a; (-1) this correction becomes large, and then the dominant contribution in this region is: 

lim G2 (w, = lim — r (2.44) 

A^o S(a;,z/) 

u;->-2a;(A/2) a;-*-2a;(A/2) 
whereby substituting Eqs. (2.38) and (2.39): 



3. High energy asymptotic behaviour of the scattering ampHtude 
3.1 The Pomeron interaction vertices 

The goal of this section is to sum over all enhanced diagrams, using a method based on the example of the previous 
section. The aim of our technique is to show that the more general diagrams for the Pomeron self-mass S (cj, z^) 
shown in Fig. 6, are equivalent to the diagram of Fig. 7 after replacing the Pomeron 1 — >■ 2 vertex with the 1 ^ n 
vertex. From a field theory perspective, when one of the diagrams in Fig. 6 is cut, a factor of 1/ (a; — J2i^{^i)) 
included in the expression for S (w, v), where the sum is over all the Pomerons in the cut, with BFKL kernel uj (z/j). 
In this approach, each cut brings an additional pole in the w-plane. E (cu, u) can be transformed to Y representation 
by the inverse Mellin transform: 



J a- 



E(y,i/) = ^ / due'^^^iuj,!^) (3.1) 



lim G2 (co, u) = lim 

A^O A-i-O 
£j-)-2w(A/2) uj^2uj{\/2) 



327r2a2 A + A+ A- A 



(A + A+) (A - A+) 



327r2a2A4 

Since A+ is small in the limit that w — >■ 2a; (A/2) (see Eq. (2.37)) then Eq. (2.45) simplifies to the following asymptotic formula: 



(2.45) 



lim G2 {cj, v) 

A— >-0 
tj-5>2tj(A/2) 



lim 

A-i>0 
a;-i.2i^(A/2) 



327r2a2 A"* 



A 

w — 2a; I — 
2 



(2.46) 



Eq. (2.46) is the dominant part of the propagator of the dressed Pomeron. The behaviour of the dressed Pomeron propagator with 
energy can be seen by transforming to Y representation using the following inverse Mellin transform: 



^dressed ^y^^^ 



/■a+ioo 
J a—ioo 



— e &(a;,.) 



(2.47) 



Eq. (2.47) leads to the amplitude of the exchange of one dressed Pomeron, that grows with energy according to the following behaviour: 



^dressed (y^^) oc 



y3/2 



2i^(A/2)y 



(2.48) 



We believe that we have learned two lessons from this re-summation. The first one is that the enhanced diagrams change the asymptotic 
behaviour of the scattering amThroughout the paper the notation lim^i-fa /(a;) is used to denote the asymptotic behaviour of f(x) at 
a; — > a. It is not meant in the conventional way as used in analysis 



Using Eq. (3.1), the contour of the w-integral can be closed over each pole that stems from 1/ (cj — w(^'i) )• 
The residue from each pole wih lead to the expression for S (y, i/) oc exp (^^ a;(fi)y) (where Y is the energy 
variable for dipole - dipole scattering.) The largest contribution to S (Y, u) stems from the pole 1/ {oj — Y17=i ^ i^i) ) 
in the w-plane, where n is the maximum number of Pomerons, that can be cut in the diagram. The residue of this 
pole yields the contribution to E {Y, v) of the order: 



where = 4(5^ In 2 is the leading order contribution to the intercept of the BFKL Pomeron, (see the expansion 
of Eq. (A. 4) and the surrounding discussion). Based on this observation, we propose the following method for 
calculating the Pomeron self-mass S(a;,i^), for the general diagrams of Fig. 6. Consider the diagram, where the 
maximum number of Pomerons in a cut is n. Then E„ (w, v) is proportional to: 



where each pole 1/ {uj — u (fj) ) stems from a different cut in the diagram. The term 1/ ( a; — Ym=i ^ i^i) ) 
comes from the cut, that cuts the maximum number (n) Pomerons. To transform to Y representation, Eq. (3.3) 
should be substituted into Eq. (3.1). In our approach, we close the w-contour around the pole 1/ { u — X"=i ^ (^i) ) 
(the maximum Pomeron cut). Then the solution is equal to the residue of this pole, i.e. we replace u = 
Sr=i i^i) everywhere in the integrand. The remaining poles are absorbed in the expression for the 1 — >■ n 
Pomeron vertex, (which we will derive below). In this way the diagrams of Fig. 6 are equivalent to the diagram 
shown in Fig. 7. In this approach, S„ can be calculated according to the following formula (which is shown 
graphically in Fig. 7): 



where {ui} denotes 2/1,1^2, .. . Un- This method of calculation, is directly related to the Mueller-Patel-Salam- 
lancu approximation, for calculating the main contribution to the scattering amplitude due to the exchange of 
BFKL Pomerons [20]. Strictly speaking, Eq. (3.4) is the t-channel unitarity constraint in (cj, z/)-representation. 

The diagrams for the 1 ^ n Pomeron vertices are shown in Fig. 8. The simplest diagram for the 1 — > 2 vertex 
shown in Fig. 8 a), was calculated in the previous section in detail. It is useful to illustrate the main steps of this 
calculation, since this approach can be easily generalized to the calculation of the 1 — > n vertex, for arbitrary n. 
We draw attention to the formula for the 1 — )■ 2 vertex given in Eq. (2.5). Recall that in this formula, u is the 
scaling dimension of the Parent Pomeron, and z^i,f2 are the scaling dimensions of the two daughter Pomerons, 
that are produced at the vertex. The dominant contribution to the 1^2 vertex stems from the singular region 
1/2 — 1/ = 1^1 + 1/2- Closing the integration contour around this pole, leads to the conservation relation: 




(3.2) 




(3.3) 




(3.4) 




Figure 6: Examples of the diagrams for tlie Pomeron Figure 7: The graphical representation for the formula 

self-energy S. Dashed lines denote the cut with the max- of Eq. (3.4), for the Pomeron self-energy I]„. 
imal number of Pomerons. Wavy lines are used for BFKL 
Pomerons. 




Figure 8: The diagrams for the multi-Pomeron vertices. Fig. a) shows the simplest 1^-2 vertex. Fig. b) and Fig. c) show 
the 1 — > 3 vertices, and Figs, d) - i) show the 1^4 vertices. The dashed lines show the cross sections with different energy 
(lu) propagators. Wavy lines describe the BFKL Pomerons. 

i/2 - V = ui+ 1/2 or Ai2 = Ai -|- A2 (3.5) 

where A12 = i/2 — z/ and Aj = Ui (i = 1,2). The values of i^i and 1^2 are small, since this leads to the dominant 
y-dependence of the simple loop of Fig. 2, proportional to exp {2ljJqY)K Hence from Eq. (3.5) we can conclude 

"The dominant contribution stems from small ui and U2, which are of the order ln{x^) / y^uj "(0) Y <C 1, where x'^ is the variable. 



that V i/2, (or in other words i/2 — v = A12 — >■ 0) at the 1^2 vertex. 



Now we generaUze to the notation f j (z = 1, 2, . . . , n), where ui denotes the scahng dimension of the n daughter 
Pomerons, produced from the 1 ^ n vertex. The values of the fj's are small, which leads to the y-dependence of 
the diagrams proportional to exp(najoy). At the simplest 1 — >■ 2 vertex (see Fig. 8 a) ) where z/12 i^i,t'2) the 
conservation law i/2 — V12 = vi + V2 holds. Hence for z^i, V2 small, then V12 = i/2. At the 1^3 vertex shown for 
example in Fig. 8 b) 1^123 1^12,1^3, and at the 1 — > 4 vertex shown for example in Fig. 8 h), 1^1234 z^i2,f34.In 
general this leads to the conservation rule at large n: 



yi2...n = i{n-l) /2 + ^\i Ai < 1 (3.6) 

1=1 

where the scaling dimensions of the produced Pomerons, are denoted by Aj, where \i<^\. Here n is the integer 
number which counts the number of Pomerons produced in the tree decay, that started with one Pomeron. Within 

the general 1 — t- n vertex diagram, i'i2...n is close to i (n — 1) /2. Note that Fig. 8 a) contains one vertex, Fig. 8 b) 
contains 2 vertices and Fig. 8 d) contains 3 vertices, such that one can generalize to the 1 ^ n vertex that contains 
ra — 1 sets of 1 — 2 vertices. The expression for the 1^3 vertex diagrams, as illustrated in Fig. 8 b) and Fig. 8 c) 
can be written as follows: 



r(z^l23|z^l,Z^2,i^3) = (3.7) 



F (Z^123 I Z^12,Z^3) 5(^^12) F (z/12 I 1^1,1^2) _^ F (1/123 | 1^23, l^l) 5 (^^23) F ( Z/23 1 Z/2 , /^3 ) 



U) — 0j{vi2) - w(f3) W — 0J{V23) - Oj{l'i) 



In Eq. (4.1), oj = Y^^^i oj (vi) because we are taking the residue of the a;-plane pole 1/ (^lo — Yl^=i ^ (^i)) 5 which 
comes from the cut in the diagrams of Fig. 8 (b) and Fig. 8 (c), that cuts all 3 Pomerons. Eq. (4.1) can be be 
written as 



F (z/123 I Z^l, 1^2, 1^3) = r (z/123 I 1^12, J^s) f {l'l2 | Z^l, 1^2) + F (1/123 | 1^23, Z^l) f (1^23 I Z^2, l^d) (3.8) 

where the following definition was introduced: 

f,/ I X 5(i^l2)F(z/i2 |z/l,Z/2) .„„^ 

F (1/12 1/1,1/2) = , . . — 7 r (3.9) 

LO [Vl) + UJ (f2) - U (z/12) 



built from the size of the dipole (see ref.[14] for the formula for x^). 



Figure 9: The graphic form of the equation for f , given in Eq. (3.10) 



with a similar definition for f (i/'23 \i'2,'^3)- The equation for the vertex where 1 Pomeron — >■ n Pomerons is 
shown in Fig. 9. The equation for the vertex takes the following form: 



n-l 



t{n,X\{Xi = iyi}) = 5^f(n-l,A|{Ai},A,) f(2,Aj|Ai,Afe) (3.10) 



where 



f(2,A,,|A„A,) = V^^ai^'iai ^^-^^^ 

This equation tells us, that the emission of one extra Pomeron (shown by a zigzag line in Fig. 10), can be 
reduced to the emission of one extra Pomeron from the produced (n — 1) daughter Pomerons. Diagrams where 
the extra Pomeron is emitted elsewhere in the decay tree, cancel. Indeed, the vertex for one Pomeron emission 
r (2, Ali^j, ffe) can be re-written in the form 

r(2,A,fe|Ai,Ait) = f (2,A,fe|Ai,Afc) { G-i(a;,sW) - G-i(a;,s(2)) } (3.12) 
where, using the example of Fig. 10 - B: 

= a; - sW = w - a; (A012) - w (A3) (3.13) 

= a;-S(2) = a; - a; (Aq) - a; (A12) - w (A3) (3.14) 

Using this Ward identity we can show, that thanks to the cancellations of the diagrams, (as shown by the 
example of Fig. 10), the emission of the extra Pomeron occurs only from the produced (n — 1) daughter Pomerons 
in the final state. Diagrams where the extra Pomeron is produced from intermediate Pomerons in the decay tree, 
cancel. In Fig. 10 we show the use of the Ward identity when calculating the 1^4 vertex, in terms of the 
1 — >■ 3 vertex. The blob in Fig. 10 is used to denote the product of two vertices r(l 2) r(l — >■ 2) without any 






Figure 10: The illustration of the Ward identity of Eq. (3.12). The extra emitted Pomeron is denoted by the zigzag Une. 



G{uj, Vik) between them. After summing all of the diagrams in Fig. 10-A,Fig. 10-B and Fig. 10-C, one can see that 
f(4,A|z^i,i/2,i^3,'^4) = Zli=ir(3,A|Aifc,z/;,^'j)f(2,Aife|z^i,zvfc) where ii^k + l + j. 

After this general outline of the calculation using Eq. (3.12), we calculate the example of r(l — )■ 4) in more 
detail. Each diagram of Fig. 10 can be written as a product of three F (2, Aijfc|Aj, Afe) terms. Since F (2, Aife|Ai, Afe) oc 
Ajjfc/(AjAfe) (see Eq. (3.19) below), one can see that this product turns out to be the same for each of the diagrams 
in Fig. 10. Having this in mind, and using Eq. (3.12) we have for the diagram of Fig. 10-B 

Fig. 10-B ^ f (4,A|Ao,Ai,A2,A3) = f (2, A|Aoi2, A3) f (A012IA12, Ao) f (A12IA1, A2) 
X G (c, { G-^ (c, S«) - G-^(a;,s(2)) 



nf4G(a;,s(2)) G-i(c.,S«) - l} 



(3.15) 



Here 11?=! is used to denote the product of corresponding f terms, since it does not depend on the diagram. 
The first term in Eq. (3.15) describes the emission from the produced Pomerons and has the same structure as the 
diagram of Fig. 10-A at a; = X^^=o^ l'^*)- Indeed, f for this diagram has the form 



Fig. 10 - A 



^ F(4,A|Ao,Ai,A2,A3) = F (2, A|Aoi2, A3) F (A012IA02, Ai) F (A02IA0, A2) 

X G (a;, S(2)) ( - a; (A012) + ^ (Ai) + a; (A02) ) t^^t 7^^^^"' [^^ ^^'^) 

V /V / a; — a; (A012) — w (A3j V / 

3 3 

U)=X;<=o'*'('^i) 



nfiG(a;,SP)) (a;,S«) 



(3.16) 



j=i 



The difference between Eq. (3.16) and the first term in Eq. (3.15) is in the expression for S^^^ which is equal to 
= u (Ai) + w (A02) + w (A3) instead of Eq. (3.14). 



The second term in Eq. (3.15) cancels with the diagram of Fig. 10-C. This diagram is actually equal to sum of 
the two terms shown in Fig. 8-h and Fig. 8-i. This sum is equal to the following expression. 



1 

Fig. 10-c ^ f(4,A|Ao,Ai,A2,A3) 



(a; (Ao) + u (A3) - uj (A03)) 
+ (a;(Ai) + a;(A2) -a;(Ai2)) 



u) - uj (A12) - a; (A03) 

u;(Ao) +a;(A3) -^(Apa) 
oj - uj (A03) - oj (Ai) - uj (A2 
oj (Ai) + oj (A2) — oj (A12) 



uj-u) (A12) - u (Ao) - a; (A3) 

3 

J{ti (3.17) 



Summarizing we see that Eq. (3.12) leads to the cancellation of the emission of the Pomeron from the internal lines 
in the diagram. Such cancellations are analogous to the cancellation in gauge theories since Eq. (3.12) is similar to 
the Ward identity in these theories. 

As discussed above, at the 1^2 vertex, the conservation rule i/2 — Vj = Ui + Vk {^ik = Aj + Afe) (see Eq. (3.5) 
and the surrounding discussion) . We are interested in the region where f j and Uk are small, since this leads to the 
y-dependencc of the diagrams proportional to exp {nojQY), which is the dominant contribution. Thus the relevant 
region is i/2 — z/j — > 0. Substituting for F {vj \i'i-,i'k ) ttie formula of Eq. (2.5), then Eq. (3.11) becomes: 



^-'^planar (^i \^i^^k) ^2 -'^nonplanar (^i \^i^^k )^ 



r{2,i^j\i^i,i^k) = ^ g{i^j) ^ r-. r'h (3-18) 



Now to calculate the residue of f (2, Vj \ Vi, v^) in our region of interest, namely i/2 — Uj — i/i — ^ when 
t^iji^k 0, simply substitute into Eq. (3.18) the asymptotic formulae for Fpiaji^r and FuQ^planar derived in Eqs. 
(A. 22) and (A. 26). Note that in this region, as i/2 — Vj ^ Vi + ^ ^ tbc denominator of Eq. (3.18) tends to 
oj{vj) — uj{ui) — uj{vk) — ^ Cis/ {\ + ivj) (see definition of Eq. (A.l) ) and g {uj) — ^ — (5 + ii^j) ^ /4 (see definition 
Eq. (2.1)). With this in mind, overall Eq. (3.18) in this region reads: 



i + iiyj) 



lim F (2, Uj I Ui, Uk) = lim 

= lim bi!^l±^ (3.19) 



(3.20) 



where the constant a is defined in Eq. (A. 28). Now Eq. (3.10) is the equation for the BFKL Pomeron fan 
diagram. Fan diagrams can be summed using the generating functional technique that has been developed in {Y, Q 
representation** (see ref. [32] for details). For the generating functional, we can write the linear equation in terms 
of functional derivatives, which reflects the fact that the Pomeron can decay into two Pomerons. In the dipole 
model, this decay can be written as the decay of one dipole to two dipoles. Eq. (3.10) simplifies this functional 
equation to a recursive formula. Two simplification rules are essential for our approach: (i) the most singular part 
of the triple Pomeron vertex has a much simpler form in v representation, than the BFKL kernel in coordinate 
representation, and (ii) the loop correction to the vertex can be neglected at high energies (see ref. [20] for a full 
explanation). Eq. (3.10) can be viewed as an equation in time (rapidity). Indeed, Eq. (3.10) states that the process 
of n-Pomeron production can be considered to be the production of n — 1 Pomerons at time t, and the later decay 
of one of the produced Pomerons into two, at time t + 6t {6t <C 1), as shown in Fig. 9. Using Eq. (3.19) we can 
rewrite Eq. (3.10) in the following form: 




(3.21) 



where on the LHS of Eq. (3.21) the notation {fj} = ui, f2, . . . I'n, whereas on the RHS the notation {ui} = 
ui,i>2, ■■■I'n', ^j- In this approach Eq. (3.21) yields the following solution: 



" 7 

V{n,\\{vi}) = \{n-l)\Y[^{vi) with $ (z/^) = - (3.22) 

It is easy to check that this solution satisfies both the recursive equation (see Eq. (3.21)) and the initial 
condition of Eq. (3.19). 



3.2 Green function of the resulting BFKL Pomeron 

Using the formula of Eq. (3.22) for the vertices, we can calculate E„ (w, v) from Eq. (3.4) as: 



n — 1)\ 1 / -TT dui 



The explanation behind the factor in front in Eq. (3.23), is as follows. The vertices that enter into Eq. (3.23), 
are equal to T after multiplying by a factor of (w {v) — u (vi) — oj (^2)) /fl'(z^) (see the definition of Eq. (3.11)). In 



**C is conjugate variable to v (see more details below) 



the region where A = z/2 — v ^ Q, this factor reduces to the asymptote a^/ {Xgiu)) for large n. The factor of 
(— l)"'(n — l)!/n that appears in Eq. (3.23) has the following meaning: 



(-l)"(n - ly./n = { (-1) for each Pomeron loop} x (n - 1)!^ { from Eq. (3.22) } 

X — { from the identity of the Pomerons} x (—1) { from definition of S} (3.24) 

We do not need to integrate over the entire phase space in Vi, due to the identity of Pomerons. It is enough 
to integrate within the region I'n > ^n-i > ■ ■ ■ > • • • > I'l- As one can see, in this region all of the Pomerons 
have different z^'s, and can be considered to be different particles. This region covers the 1/n! part of the entire 
phase space in v. Further summation turns out to be simpler in Y and ln(^ = Inxx* representation, where Y is the 
rapidity of the dipole-dipole scattering, while 



where b is the impact parameter of the dipole-dipole scattering, and f and R are the sizes of the two dipoles. 
Using these variables, E„ {Y, () can be calculated using the following transform: 



'>a+ioo fta+oo 



^n{Y,0 = I l^e-^/ ^e-i-^S„(c.,.) (3.26) 

Jia'- 



2m / , 27r 

' a— too J xa — oo 



First we switch to Y representation using the following inverse Mellin transform: 



J a—ioo ' J i \ i / 

Then we switch to C, representation using the following approach: 



/»ia'+oo 
J ia' —oo 
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1 ^ 3 
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1 ^ 6 
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40 



Table 1: The number diagram for tlie values of z/„ for diffrent vetices F. 



As we have discussed generally speaking v = — A^^. However, the value of turns out to be different from 
(n — 1)! in general, for various different 1 — )■ n vertex diagrams. In Table 1 we give the examples for 1 — >■ n vertex 
diagrams, up to n = 6. 

Unfortunately, we have not derived the general rules for how to calculate the value of Vn for the 1 ^ n vertex. 
We consider two models for such numbers: (i) each \ ^ n vertex has = 0, and (ii) each 1 — ^ n vertex has 
i^n = i(n — l)/2. The first model gives the sum of the leading twist contribution (^'^ ^ 1, while the second model 
sums over all high twists C~"'^^- We believe that considering these two models for finding the value of f„, 

provides the largest possible contributions. This belief is based on the following simple examples. As can be seen 
from Eq. (3.28), we are summing an asymptotic series of the form: 



oo 

^(-irC„L" (3.29) 

n=0 

where L is a large parameter. Our first model means that for the leading twist contribution, we choose Un = ^ 
in all n!-diagrams. This leads to Cn oc n\. In the exact approach, the number of diagrams with i/„ = is less 
than nl. However, the largest sum corresponds to C„ = n\. One can see this by setting C„ = 1 and C„ = 1/n! in 
Eq. (3.29). The same occurs in the second model, which we believe leads to the maximal sum of the highest twist 
contributions. 



3.3 Summing high twists 

Recall that u = i(n — l)/2 + ARii(n — l)/2at large n^^, such that after switching to the integration variable 
A, then Eq. (3.28) can be written as: 



Jia'—oo J 1=1 \ i / 

i=l 



t^We recall that A = \n(^/asY < 1 
"We recall that A = In C/oisY < 1 



where we took the integral over A using the (5-function. Using Eq. (A. 4) we can solve the integral over i/j 
explicitly, using the method of steepest descents. In this approach Eq. (3.30) simplifies to: 



= ^ (-l)''"^-^^^^^^r(n) L" (3.32) 
where L = 6%/^ exp fa;(0) F - -J^") (3.33) 



Using the integral representation for the Euler-Gamma function (see formula 8.310(1) of ref. [33]): 



r(n + i) = r ee-^dt- 

Jo 



(3.34) 



and since we are summing from n = 2, (since the first 1 — )■ n vertex in the sum is the 1 — >■ 2 vertex), we obtain 
the following result for S {Y, Q: 



n=2 J n=2 

«^ /7 /"°° dt _j ( Lt + 2{Ltf - ln(l + Lt) - 2Lt ln(l + Lt) - {Ltf ln(l + Lt) 



'0 

"00 



_ai r- I dT _2./j^r + 2r2-ln(l + r)-2rin(l + T)-T2ln(l + T) 
= f VCS(L) (3.35) 



where T = Lt. For large L we find that 



fe-^/^(2-lnT + o(i^)) « f (- InT) = - (lnL)V2 (3.36) 

Inserting Eq. (3.36) into Eq. (3.35), using the definition for L given in Eq. (3.33), we see that at large Y and 
in terms of In (, then E {Y, Q tends to the following simple formula: 

m.O = f VC Ur - ^ - ilny -i.n( + ^n(„^J^))' (3.37) 



Using the following formula to transform to (w, ) representation: 



/*oo /*oo /*oo noo 

Jo ^0 ./-oo ^0 



(3.38) 



Note that Eq. (3.38) is the double Mellin transform from {Y, () to (w, v) representation. This corresponds to 
the inverse-Mellin transform of Eq. (3.26) which was used to transform from {oj, v) to (Y, Q representation. Then 
inserting Eq. (3.37) into Eq. (3.38) for large F, leads to the following equation for S (a;, v): 

S (., .) $ +"°"^7 + ^("')) (3.39) 

where the last term ©(a;^), labels small terms proportional to w^. In reality, due to the a;^ term in the 
denominator of Eq. (3.39), then overall this O (w^) term, leads to terms which are singular in powers of just l/w, 
and hence smaller than the singular terms by powers of X/uP and \/u^ that come before. The Green function of 
the dressed Pomeron reads: 

G{ijj,v) = , } ^, = 5 7 5 (3.40) 



For small uj and we can neglect the contribution from Gq^ {uj, u) = uj — uj {v). Therefore, in this limit the 
Green function of the dressed Pomeron is given by the following formula: 



8 a;3 /I 

^ ^ ^0 V2 / (3 4^^ 

where in Eq. (3.42) the second term in the denominator of Eq. (3.41) is neglected in the limit that a; — >■ 0. One 
can see that Eq. (3.42) leads to the function that has no singularity at — > and therefore, the corresponding 
imaginary part of the scattering amplitude vanishes. To gain a better understanding, keeping the second term in 
the denominator of Eq. (3.41) yields: 



Passing to {Y, (^)-representation we see that the asymptotic behaviour of the Pomeron Green function is 



■»ia'+oo ra+ioo 
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Closing the contour of integration on the real negative w-axis, we have: 



2 ,,,3 



at Lo, 







v2 /, I? 
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2 ^ ainC 
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2 ^ ainC 



J(lnC) / 
Jo 

4! 



do; e 



(3.45) 



Prom Eq. (3.45) its clear that the Green function leads to the cross section that decreases as and the 

character of this asymptotic behaviour depends on (. The factor that depends on ^, in front of Eq. (3.45), indicates 
that the main contribution turns out to be a leading twist contribution. 



This decrease at ultra high energies is the most salient result of this paper. It is well known that the Pomeron 
calculus in zero transverse dimensions, leads to the Green function of the Pomeron that decreases with energy, (see 
refs. [21, 22, 23, 24, 25]). However, in our case the decrease of the cross section is only logarithmic, whereas in the 
Pomeron calculus in zero transverse dimensions, the cross section falls exponentially at large Y. Therefore, we claim 
that the BFKL evolution of the sizes of the interacting dipoles, do not lead to a new qualitative effect. However, 
it is interesting that the character of the asymptotic behaviour, crucially depends on the size of the interacting 
dipole. 



3.4 Summing the leading twist contribution 

In the previous subsection we demonstrated that summing the contribution of high twists, leads to the solution 
which appears to be the leading twist contribution. Therefore, it seems reasonably possible to derive the leading 
twist contribution. We will do this assuming that i^n = for every odd n (see Table I). Recalling that g{u) 16i^^ 
as ^ 0, and using the formula = f„ — A then we can re-write Eq. (3.28) in the form 



/•m'+oo 
J ia'- 



la — oo 



Pia +00 r n / n \ 

- ^li-m-'^ ^ |_l,e"'"c/n^e"<".)., ,3.40) 

«/m — 00 J i \ i / 

Using the following representation for the 5- function, namely, 

(5^A-^z^i^ = ^f (i/xe^'^(^-^?=i'^0 (3.47) 



27r , 

i=l / J —00 



and integrating over A by closing the integration contour around the double pole at A = 0, we obtain 
S.(y,C) = §(-lr&^"^^^^^/^(^l^C + ^/^)^/^e--''^ + -(^')^ 

J —00 i=l J 

(3.48) 

Integrating over vi using the method of steepest descents, we have: 

T.AY,Q = |(-ir6-(!i^±y^^d^(il„<+iM)(y^exp(.„l-- (3.49) 

/oo 
dxexp (— ax^) = \/i^/a yields: 
-00 

S„(y,C) = -flnC(-l)"^ (3.50) 



16 n V vrn 



L = tS/^e-W^' (3.51) 



Using Eq. (3.34) we can re-write Eq. (3.51) in the form 



n=2 



'00 



U n=2 



-2 / /""^ 

^y^lnc/ ^e-*(Lt + Lz3/2(-Li)) (3.52) 



For properties of the Polylogarithm function see Ref.[36]. Since we assumed that Un = occurs only at odd 
n, we can extract this from the sum of Eq. (3.52), by subtracting from it the function where L — )■ —L. In this 
approach, finally we arrive at the expression: 



Oil DY 



16 V TT 



InC 2L + 



dT 



exp (-T/L) (Lz3/2(-r) - Lz3/2(r)) 



(3.53) 



16 V TT 

where T = Lt. Since the asymptotic behaviour of the Polylogarithm function Lis{—T) is known, namely 



(3.54) 



Lis{-T) ^ 



In' (T) 



r(i + s) 

Then with this in mind, the integral of Eq. (3.53) is expected to lead to the following result: 



(3.55) 



dT 



S(L) = 2L+ exp{-T/L) {Li:,/2{-T) - Lig/s (T)) 



3/2 



L»i dT (InT) 

^ Jo ^'"^^"^/''^TW^ 



2L 



2 (In L) 



5/2 



5 r(5/2) 



(3.56) 



Therefore at large F, after inserting Eq. (3.56) into Eq. (3.54), yields the following high energy behavior for 
S(F,C): 



5/2 • 



(3.57) 



Using the formula of Eq. (3.38) to switch to (a;, v) representation, by approximating the second term in 

5/2 



Eq. (3.57) at y as (^qY - ^ InY + In (^"^ 

the formula: 



^ (^^Y^ - ^ liyYu^Y'^ + O (y)^ , leads to 



E(a;,z/) = 



q2 1 



26^ 



+ 



16 I a; - Wo r (7/2) 



a;4 



5 a;, 



3/2 



2 



(Inw + ©(Const)) 



(3.58) 



where O (Const) denotes constant terms that do not depend on uj. Plugging Eq. (3.58) into the formula for 
the dressed Pomeron propagator G2 (w,z^), defined in Eq. (2.29), then one arrives at the following expression: 



G2 (w, u) 



oj — oj {u) + S (w, v) 



U) — (jj {v) + 




(3.59) 
(3.60) 



In w 



From inspection of Eqs. (3.58) and (3.59), then uj ^ uoq leads to the asymptote G2 (i^, J^) cc uj — ujq and 
therefore, this region does not contribute to the cross section due to the absence of singularities in the region 
UJ = uq. In the limit that w — > 0, Eq. (3.60) tends to the following limit (where the third term in the denominator 
of Eq. (3.60) dominates in this region): 




The first term in brackets in Eq. (3.61) leads to G2 {uJ, oc uj^, which also does not contribute to the total 
cross section, thanks to the absence of any singularities in this region. Hence it follows that the first contribution 
stems from the second term in Eq. (3.61), which leads to the contribution to G2 {oj, u) equal to: 



G2{oj,u) = ^ Ino. (3.62) 

where in passing from Eq. (3.61) to Eq. (3.62) the fact that F (7/2) = 157rV2/8 was used. Finally passing to 
(Y, C)-representation we see that the asymptotic behaviour of the Pomeron Green function is 



I '+00 /*a+ioo 



G2{Y,0 = r '"^e^'^'^i r""^e-^G2{u;,,.) 

J ia'- 



27r / . 2m 

'—00 J a— too 



Hl-C)^ (3.63) 



Hence its clear from Eq. (3.63) that the summation of the leading twist contribution, leads to the same 
qualitative result, namely that the Pomeron Green function vanishes at large Y, but only logarithmically. This 
style of decrease is steep enough to provide the final answer, without the need for further re-summation. 

It should be mentioned that in spite of the decreasing behavior of the amplitude with energy, at large values 
of b the integral over the impact parameter turns out to be divergent, indicating that the problem of the large 
6-dependence cannot not be cured by summing enhanced diagrams. Moreover, this problem needs new ideas from 
non-perturbative QCD, in order to find a solution. 



4. Conclusions 

This paper describes the technique developed to find the sum of enhanced diagrams ( Pomeron loops), in the 
dipole-dipole scattering process. In conclusion we would like to mention two main features of the result. The first 
one, that the cross section and/or the Green function of the dressed BFKL Pomeron falls down with energy. The 
second result, is that the asymptotic behaviour depends crucially on the size of the colliding dipoles, and the impact 
parameter of the collision. 

We wish also to draw the attention of our reader, to two selection rules which are essential for our approach 

to the summation of the enhanced diagrams. First, we restrict ourselves to the contribution to the triple BFKL 
Pomeron vertex (see Eq. (2.5)), that is singular in the region i/2 — ly — — 1^2 — > 0, for small i^i and 1/2, since 
this part of the vertex generates the most singular contribution, leading to a larger result than the other parts of 
the vertex. Second, we neglected the contribution from Pomeron loops to the dressed vertices. This assumption 
is equivalent to the Mueller-Patel-Salam-Iancu approximation [20], formulated in the s-channel of the reaction. 
Therefore, we sum the Pomeron loop diagrams, but using the above mentioned specific assumption about Pomeron 
vertices. It should be stressed that the Mueller-Patel-Salam-Iancu approximation, as well as our approximation, 
selects the diagrams with the most essential increase with energy, and therefore, it can be used in our approach. 

As we have discussed in the introduction, we cannot prove the BFKL Pomeron calculus, based only on the 
triple Pomeron vertices. Moreover, we personally have an argument stating, that it is necessary to introduce the 
four BFKL Pomeron vertex [18]. The fact that we obtained the total cross section that decreases with energy, 
stands as a reminder of the constant cross section obtained, only after taking into account the four Pomeron vertex, 
in 1 + 1 dimensional Pomeron calculus. However, we would like to stress that in our case, there is only a logarithmic 
decrease, which is steep enough to claim that the summation of the enhanced diagram provides the solution to the 
problem. 

As we have mentioned in the introduction, using the BFKL Pomeron calculus in the form of Eq. (1.1), we 
neglected both the vertices of transition of one Pomeron to more than two Pomerons, as well as the contribution 
of the multi -gluon states in the next-to-leading 1/iVc order (see Fig. 11). 

In Fig. 11 the wavy lines denote the dressed Pomeron Green function (G), which is proportional to or to 
depending on the model. Therefore, the contribution of the diagram of Fig. 11-a is equal to 



Y 




r - 

0-5. 



a) 



b) 



c) 



Figure 11: The vertex of one Pomeron to three Pomerons (Fig. 11-a); the contribution of the multi-Pomeron states (Fig. 11- 
b) and the first correction to the two Pomerons state (Fig. 11-c). The wavy hnes denote the Green function of the dressed 
Pomeron. 



One can see that the new vertex F (1 — > 3) led to the renormalization of the coupling of the Pomeron with the 
dipole, that is included in Eq. (1.1) in the term Se- In other words, the new vertex did not change the BFKL 
Pomeron calculus. 

The multi-gluon state, i.e. 2n-gluons in the t channel, generates a larger intercept, than the intercept that the 
corresponding number (n) Pomerons would generate. As shown in Ref. [19], the new intercept can be calculated 
by summing the Pomeron exchanges, with the new vertex of the interaction of two Pomerons (see Fig. 11-b). The 
contribution of the first diagram shown in Fig. 11-c is equal to 



Therefore, the contribution of the multi -Pomeron exchanges with the interaction between them leads to small, 
negligible contributions. 

In light of this, we conclude that the correction to the BFKL Pomeron calculus cannot change the fact that the 
total cross section falls at high energy. Of course, we have discussed only those multi-gluon states in the t - channel, 
where each pair of gluons is in a colourless state. These states have been discussed in Ref. [19]. As far as we know, 
this is the only known example of the source of the intercept, which is larger than the intercept of the exchange of 
several Pomerons. Unfortunately, we do not know about the situation with the exchange of more general n-gluon 
states, which could make the BFKL Pomeron calculus incorrect in the 1/Nc approximation. However, the examples 
of such states have not been found ( see Ref. [37]). 




^ r (1 ^ 3) / dY' {Y - Y') ^ Y'-^^ oc r (1 ^ 3) oc G (Y) 



(4.1) 




^ r (2 ^ 2) / dY' (y - y')~'" y'-i° oc r (2 ^ 2) y-^" oc (y) < g (y) 



(4.2) 



We would like to recall, that the final re-summation was applied in the case of two extreme models, for the v 
dependence of the 1 Pomeron — )■ n Pomerons transition vertices. The answer to this problem, is the next challenging 
problem which we hope to resolve in the near future. 

Nevertheless, we believe that the problem addressed in this paper, is the first to be solved in order to start a 
theoretical discussion on what could happen with the dilute-dilute system of scattering, at high energy. Observing 
the substantial difference between the Pomeron self-energy and its Green function, we conclude that the quantum 
effects due to the BFKL Pomeron interaction, change the character of the asymptotic behaviour. Since these effects 
are so strong, they have to be taken into account for dilute-dense and dense-dense scattering, at ultra high energies. 
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6. Appendix A - The BFKL kernel and the triple Pomeron vertex 



The BFKL kernel oj {v) is defined as 



^(z^) = OisX{^) \ i^Ots = (A-1) 

X{y) = + = 2^{l)-^{^ + iu^-^{^-iu^ (A.2) 

where tp (x) is the Di-Gamma function, defined as [33] : 

*W = = -^--f(^-ril) (A.3) 

In the vicinity of i/ — > 0, the BFKL kernel takes the following form 

u{y) = 4as ln2 - 14a5C(3)z^^ = a; (0) - Du"^ = ojq - Du"^ (A.4) 

The planeir vertex 

The triple Pomeron vertex T [u \ ui, 1^2) was defined above in Eq. (2.5), in terms of two contributing diagrams, 
namely the planar and non-planar diagrams shown in Fig. 3. The expression for the planar diagram given in 
Eq. (2.6), contained the function Q {y\vi,U2)., which is given by the following integral (for a detailed explanation 
and derivation of these results, see Ref. [34]): 



J \ Z01Z12Z20 \ \Z0Zi ZqZi J Zi){l- Z2){1- Zi){l- Z2) J \ J 

The integrations of Eq. (A.5) were evaluated in [34], where the following results were derived. Note that in the 
original paper of ref. [34], the symmetric property ft {vi \ 1^2, 1^3) = ^^(^"3 1 1^1,1^2) = ^(2^2 1 t's, z^i) was derived. We 
have used this symmetric property, and hence the version of the function J7 written below is based on a different 
permutation of arguments and differs from the version given in ref. [34] . 



fl{u\ui,U2) = Oi (z/ I 1/1,1/2) + ^^2 (i^ I 1^1, ^'2) + ^^3 I 2^1,1^2) 



(A.6) 
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J^3(z/|z/i,Z/2) = fl2{v\v2,Vl) 



(A.9) 



Using the following identity for the Hyper-geometric function (see Ref. [33] formula 9.100) 



2-Fi ^ a, 6 c = (1 — " ^ 2F1 (^c — a, c — b c 



(A.IO) 



then the factor 2^1 ^— ^ "ii^i, — — ii^i jl 
may be re-written, such that: 



X ] 2F1 { - — ii'2-, - — \\\x \ that appears in Eqs. (A. 7) and (A. 8) 
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(A.12) 



J]3(zy|z/l,f2) = ^2 (l' I 2^2,2^1) (A.13) 

Using the fohowing expansions for the Hyper-geometric functions (see Ref.[33] formulae 9.100 and 9.14): 
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Then Eqs. (A. 11 - A. 13) can be re-cast in the following form: 
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S(m|n) = / dxx'^-Ul-xT-^ = ^ ^J"^ ^ ^""^ (A.19) 
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^2,(y\vx,V2) = ^2 {1^11^2,1^1) (A.18) 

where B {m\n) is the well known beta-function, defined as (see Ref.[33] formula 8.380): 



In the region z/i — > and 1^2 — >■ we can simplify the expressions for the Oj. All of the contain a pole at 
i/2 — u — u\ — ii'2 = 0. In vicinity of this pole: 

B ( + iu + iui + iu2] -j ^ + {non singlular terms} (A. 20) 

Using the limit of Eq. (A.20) we can calculate the residue of the pole at i/2 — u — ui — 11/2 = of the function 
^2 ('^ I z^i, ^^2 ) = ^3 I ^^2, z^i ); for the even narrower region where ui and 1^2 0. In this approach, from 
Eq. (A.17) in this region, the residue is: 



lim (i/2 - 1/ - i/i - ) O2 (1/ 1 1^1, ^2) 
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1/1,1^2^0 2>Vil'2 

The identical result is true for Vt^ {v \ vi, 1/2). Repeating this procedure we obtain that f^i oc {vi + V2) <^ ^2- 
Thus in the region i/2 — v — vi — if 2 = when 1/1,^2 0, then the function 17 (z/ | z/i, 2/2) = Yl^=i I '^i, ^2) 
tends to — 27r/3z/if2, neglecting the contribution from fli. Inserting this into Eq. (2.6) leads to the result: 



lim (i/2 -u- VI -1/2) Tpianar i^' I 1^1,^2) 

i/2—i/—ui—iu2^0 

Vl,V2^Q 

The non-planar vertex 

The formula for the triple Pomeron vertex of Eq. (2.5) contained also the contribution of the non-planar 
diagram, given by Eq. (2.7). This expression included the function K{y\i'i,U2) , which is given in terms of the 
following integral in ref. [34]: 



lim !^ '— (A.22) 

Vl,V2^Q 



/ \ zoi\^ \ZQZiZoziJ V(l - ^0) (1 - ^1) (1 - 2=0) (1 - 2:1)/ V / 



The integrations of Eq. (A.23) were evaluated in ref. [34], where the following result was derived: 



A{i^\ui,U2) = 22''^+2i.i+2i.2-i^2 yV^ /_ (^24) 



Using the above expression it is easy to see that: 

lim {i/2-u-ui-u2) A{u\i'i,u2) = lim ^ '^^^'^ (A.25) 

z/2-i^-i/i-W2-*-0 i/2-v-i'i-iv2^0 2 1/1I/2 

Plugging this result into the definition of the non-planar diagram of Eq. (2.7), we can derive the non-planar 
diagram in the limit z/2 — u — ui — iv2 ^ Q where ui,V2 0. Using the fact that limj/2-iy-+o X (^) ~^ 1/ (5 + '^^) 
(see the definition of Eq. (A. 2)), the asymptote is: 

^2 a+iyf 

hm (i/2 -U-V^-V2) r lanar (l^ I ^^1, 1^2) = lim -TT — (A.26) 

Finally, inserting the asymptotes of Eqs. (A. 22) and (A.26) into Eq. (2.5), we can calculate r(i^| 1^1,^2) = 
{i/2 — v — vi — V2) r (y 1 1^1, 1^2) in the region z/2 — v — v\ — iv2 where v\, ^2 — >■ as: 

lim f(i/ 1 1/1,1/2) = lim ^—^ — '-^ {A.27) 

I/l,i/2^0 V\,V2^Q 

where „ = - ^ - + ^ j (A.28) 

In terms of the variables { A, cr, A } defined in Eq. (2.8), then Eq. (A.27) becomes: 

lim f(A|<,,A) = lim (,^;;;_^) (A.29) 

(7->-0 £7->-0 
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